Full Current Statistics of Incoherent "Cold Electrons" 



O 

o 

(N 



(N 



Dmitry A. Bagrets 
Institut fiir Theoretische Festkorperphysik, Universitdt Karlsruhe, 76128 Karlsruhe, Germany. 

(Dated: February 2, 2008) 

We evaluate the full current statistics (FCS) in the low dimensional ( ID and 2D ) diffusive 
conductors in the incoherent regime, eV S> -Eih = D/L^ , Sxh being the Thouless energy, ft is shown 
that Coulomb interaction substantially enhances the probability of big current fluctuations for short 
conductors with i?Th2>l/T_B, te being the energy relaxation time, leading to the exponential tails 
in the current distribution. The current fluctuations are most strong for low temperatures, provided 

^Th ~ [{(^V)'^ / Duf] for fD and i^xh ~ (eV^/ff) lug for 2D, where <; is a dimensionless conductance 
and vi is a ID density of states. The FCS in the "hot electron" regime is also discussed. 

PACS numbers: 73.23.-b, 72.70.+m, 05.40.-a, 73.50.Td 
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The influence of the Coulomb interaction onto the 
transport properties of low-dimensional diffusive systems 
has been a subject of extensive research for more than 
twenty years [l|, y] ■ Initially, the conductance only was a 
main object of study. The powerful alternative to that is 
to investigate a quantum noise yjj] , or, more generally, 
the full current statistics (FCS) • 

In the short diffusive wires with Etyi ^ eV, where 
ETh — D/L? is the Thouless energy, the shot noise 
equals to S* = 2\e\IF, F = 1/3 being the Fano fac- 
tor [a, lD| • In this case the conductor is coherent and ef- 
fectively zero-dimensional so that all effects of Coulomb 
interaction come from the external electromagnetic en- 
vironment. It modifies the conductance, noise '8| and 
generally the FCS d [13. 

Much less is known about the role of Coulomb inter- 
action onto the FCS in the quasi- ID and 2D diffusive 
systems, when i?Th <C eV. Under this condition the in- 
elastic electron-electron scattering inside the conductor is 
important. This subject has recently attracted the atten- 
tion in Ref . [IJ [ij, [13 , where the " hot electron" regime 
in the FCS was discussed. In this regime i^xh "C 1/te, 
Te being the energy relaxation time, and the electron dis- 
tribution function relaxes to the local Fermi distribution. 
This changes the Fano factor F from 1/3 to V3/4 jl4| . 
that was confirmed experimentally 13 ■ 

The microscopic theory [l5j of electron-electron inter- 
action in low-dimensional disordered conductors predicts, 
however, two different time scales, r^ and te, where 
T0 <C tb is a dephasing time (See Table I). It is usually 
believed [ig that classical phenomena described by the 
Boltzmann equation are governed by te ■ While the time 
T0 manifests itself in essentially quantum-mechanical 
phenomena. Since the FCS is a classical quantity one 
might expect it to cross over between the coherent and 
the "hot electron" regime on the scale Etyi "^ ^/te- 

In this paper we consider the FCS in the low dimen- 
sional ( d = 1, 2 ) diffusive conductors, taking into ac- 
count the Coulomb interaction. We show that the time 
Te is indeed responsible for the smooth crossover between 
the coherent and the " hot electron" limits if one consid- 



ers the noise and the 3 cumulant. However it is not the 
case for the higher order cumulants of charge transfer in 
the shot noise limit eV ^ T. Moreover, in this limit 
the smooth crossover in the FCS does not exist. The 
Coulomb interaction drastically enhances the probability 
of current fluctuations for short conductors iSxh ^ 1/tb. 
In this regime the higher order cumulants are given by 
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where (n) ^ 1 is the average number of electrons trans- 
fered, 5 ^ 1 is a dimensionless conductance, to* ~ 
max{_ETh,T, e*} and the scale e* reads 



*(V)^I {eV)yg^ETh, ID 

\ eVeM-gETh/eV}, 2D. 
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Eq.Q shows that each {k + l)-th cumulant of charge 
transfer is parametrically enhanced versus the fc-th one 
by the large factor eV/ui* ^ 1. It also follows from Eq.QJ 
that the higher cumulants grow with increasing the volt- 
age at ETh > I/t* and decay at E^h < I/t*, where the 
new time scale t* (eV, T) is parametrically smaller than 
Te, t* ^ Te. (See Table I). The current fluctuations 
are most strong, provided T ^ i?Th ~ l/''>(^)- There- 
fore at the strongly non-equilibrium situation the time t^ 
rather than te governs the crossover in the FCS between 
the coherent and the "hot electron" limits. 

Model and the effective action. We consider a quasi- 
one-dimensional (ID) diffusive wire of a length L and a 
quasi- two-dimensional (2D) film of a size L x L, with di- 
mensionless conductance g ^ 1 and diffusion coefhcicnt 



TABLE I: The electron scattering times for low-dimensional 
(ID and 2D) diffusive conductors, E = maxjT, eV}. At T < 

T-\V),wegetT'=r^{V). 
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D. They are attached to two reservoirs with negUgi- 
ble external impedance which are kept at voltages ±.V/2. 
The current flows along z direction. We assume the in- 
coherent regime, max{ey, T} ^ i^Thj and disregard the 
possible electron-phonon scattering, so that L <^ Lo-ph- 

Our goal is to evaluate the cumulant generating func- 
tion (CGF) S{x)- The Fourier transform of exp(— S*) with 
respect to the " counting filed" y gives the current prob- 
ability distribution P{I) (See |5|). The derivatives of S 
give the average value of current, shot-noise and higher 
order moments {{n )) of charge transfer during the obser- 
vation time to. 

To evaluate the CGF, taking into account the Coulomb 
interaction, we have used the Keldysh technique and cm- 
ployed the low-energy field theory of the diffusive trans- 
port ^17;] with the action 
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Here A — diag(Ai(i, r), A2(i,r)) is the 2x2 matrix in 
Keldysh space, where Ai 2 stand for fluctuating vector 
potentials in the conductor. Hereafter we use a longitu- 
dinal photon field, curl A = 0, thus neglecting the rela- 
tivistic effects. 

The matrix (3(r,ti,i2) accounts for the electron de- 
grees of freedom and obeys the semi-classical con- 
strain Q(r) o Q(r) = 5{ti — t-i). The action F de- 
pends on X via the boundary conditions imposed on 
the field Q at the boundaries with the left(L) and 
right (R) reservoirs [l^: Q|r-B = ^R ^^"^ ^L-l ^ 
Gl(x) = exp(ixT3/2)GLexp(-ixf3/2). Here Gl.r arc 
the Keldysh Green functions in the leads. 

With action JSJ the CGF should be evaluated as a 
path integral over all possible realization A1.2 and Q. 
We proceed along the lines of Ref. 17] and employ the 
parameterization Q = e^^Ge-**^, WG + GW = 0. Here 
field W accounts for the rapid fluctuations of Q with 
typical frequencies lo ^ i^xh and momenta q S> 1/^, 
while (?(£, r) is the stationary Green function varying in 
space on the scale ~ L. First wc integrate out the field 
W in the Gaussian approximation to obtain the non- 
linear action F{x, G, A) of the screened electromagnetic 
fluctuations. We keep only quadratic terms to F that is 
equivalent to the random phase approximation (RPA). 
At the second step one can integrate the photon field A 
and arrive to the effective action i^off[x,G]. Then the 
saddle point approximation, (5i^off[x, G]/(5G = yields 
the kinetic equation for G(e,r). 

For the rest we restrict consideration to the universal 
limit of a short screening radius t~ 
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FIG. 1: The sketcli of voltage dependence of the 6"* cumulant 
of charge transfer. Left plot - ID, conductance g — 100, right 
plot - 2D, g = 10^. The temperature changes from up to 
down, T/ETh = 1, 2, 4, 8. 



y^eV/D. In this limit we get the answer 
FMX,G] = 'jgL^-'' I d'r I ^Tr (vG.ir))' + (4) 
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where V^ is 2 x 2 matrix operator in Keldysh space, cor- 
responding to the non-equilibrium diffuson propagator: 
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with 7° = i, 7^ = fa. To derive the action Q we have 
used a local approximation, i.e. we neglected gradient 
corrections proportional to (VG ^ 1/i) ^ VVF. 

Minimizing the action Fpff under constrain G(e, r)^ = 1 
one can obtain the non-linear kinetic equation for G(e, r) 
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is the matrix collision integral. This kinetic equation 
should be supplemented by the x-dependent boundary 
conditions at the interfaces with the leads. The action 10] 
EJ is one of the central results of the paper. In the 
absence of "counting" field our matrix kinetic equation 
reduces to the standard kinetic equation for the distri- 
bution function with a singular kernel K{uj) oc w'^/^-s 
(See [la, lla Il3)- We also note that only real inelas- 
tic processes with energy transfer lo < maxjeV, T} con- 
tribute to the action (0J. In the rest of the paper we 
consider the shot-noise limit eV ^ T only and proceed 
with 

Incoherent "cold electron" regime, -Exh 2> ^/te- In 
this regime the collision term in kinetic equation is small 



and one can try to find the Green function perturba- 
tively around the coherent solution which obeys the Us- 
adel equation V^ (g°{z)V^G°{z)\ = 0. In the first or- 
der in l/{EThTE) the CGF can be found by substituting 
G'^ in the action (0J. The main contribution conies from 
frequencies T < lu < eV. After some algebra we obtain 
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Here 9^ ~ ln(u + \/ii? 
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n(x,^) = -AL^{z)R^{z)e^^[l-L^{z)-R^{z)- 

[zL^{z) + {l-z)R^{z)\{e^^-l)] (9) 

where A'lj = {eV/\u!\ — l), L^{z) = sinh(l — z) 6';^/sinh0^, 
R^iz) = sinhz 0^/ sinhfl^ and to* = max{_ETh,T}- The 
presence of Thouless energy in the low frequency cut-off 
uj* is due to the fact that the lowest allowed momenta q in 
the diffusion propagator equals to 1/L, while T takes into 
account the smearing of a step in the Fermi distribution. 
The result with the above defined cut-off u* ceased 
to be valid at sufficiently high voltages. Indeed, substi- 
tuting a zero order distribution function /o(e) = (1 — 
z)fF{(- — e.V/2) + zfpif- + eV/2) to the colHsion integral, 
one estimates the I*'* order correction 
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if e± = |e ± eV/2\ < eV and e± > max{£;Th,r}. By 
virtue of Pauli principle this correction may not exceed 
unity, ^/(i) ^ 1, which is true for e± > e* only, where the 
scale e* is given by Eq. ^. Therefore a simple perturba- 
tion theory is valid provided e* < max{_ETh)2^}- Resolv- 
ing this condition we obtain that it is the case of relatively 
short conductors, or small voltages, E^i^ > 1/r*. 

The time r* decreases with increasing the voltage and 
the l'^* order perturbation theory finally breaks down at 
Et\i < 1/t*. However, in this situation we can obtain the 
result up to the factor of order of unity using the cut-off 
oj* ~ e* in Eq. ©. The point is that the role of higher 
orders terms in the perturbation series is to smear the 
step in the distribution function /(e) around e = ±eV/2 
on the scale e* , while at e S> e* the l**' order perturbation 
theory is still applicable. Therefore the overall effect is 
similar to the increase of the temperature in the system. 

The result Q H EJ with the cut-off to* = 
max{£'Th, T, e*} enables to evaluate all irreducible cumu- 
lants {{n'')) = i^{d^ /dx^)S{x) of a number of electrons 
transfered. The expansion of n(x; z) in x starts from x^- 




FIG. 2; The log of probability to measure the big non- 
equilibrium current fluctuations {eV S> T). Curve (1) - coher- 
ent regime, (2) - incoherent "cold electron" regime, 7 = 0.2 
(3) - "hot electron" regime. 



Thus there is no correction to the current on the clas- 
sical level. The interaction correction to the noise and 
to the S'' cumulant is small by the parameter 1/{E't\^te) 
and it is dominated by inelastic collisions with the energy 
transfer u ~ el^. On the contrary, the leading contribu- 
tion to the higher order cumulants is due to Coulomb 
interaction and it is dominated by quasi-elastic collisions 
with small energy transfers (jj* < lj ^ eV . Up to the 
numerical constant the result is given by Eq. (Q. The 
sketch of the voltage dependence for the 6"^ cumulant at 
different temperatures is shown in Fig.l. The cumulants 
grow with voltage in the range Et\^ > 1/t* and decay 
if ETh < 1/t*. The enhancement of higher order cumu- 
lants are most strong at small temperatures T ^ Et\i . In 
this case their maximum occurs at eV/ETh ^ g for ID 
and at eV/E^i^ ^ g/lng for 2D. 

With result Q |HJl we can also explore the current 
probability distribution P{I) = {2tt)~^ J_ exp{— S'(x) + 
i{Ito/e)x} in the long time limit (Ito/e) ^ 1. The 
action S{x) has two branch points x — i*7j where 
7 ~ {uj*/eVy^^ ^ 1. The points ±i"f give two 
threshold currents, I^ = {e/to)dS/dx\ ^^,^ , which read 
(/± -(/))/(/) =±7/3. 

Provided / < I < I^ the probability can be evalu- 
ated with a saddle point method. Due to the smallness of 
parameter I/Ej^i^te we found that P{I) only slightly de- 
viates from the probability Pq (/) of current fluctuations 
in the non- interacting limit. For larger current fluctua- 
tions the potential ^(x) = ~S{x) + *(-^^o/e)x does not 
possess the saddle point any more and one should use 
the contour Co of a zero phase, Imr2(x)| p^ = for the 
asymptotic analysis of the integral P{I). This contour 
is pinned by the branch point x — i*7i that yields the 
exponential tails in the current probability distribution 



P(/)«exp{-5(±7)-7|/|foM, I<rorI>I+ (11) 

The results for the probability distribution are displayed 
in Fig. 2. The Coulomb interaction does not affect the 



Gaussian fluctuations. However the tails of P{I) drasti- 
cally differ from those in the absence of interaction. 

The FCS of this type can be understood as the statis- 
tics of the random current of electron-hole pairs which 
are excited by the low frequency fluctuations of the elec- 
tromagnetic field, produced by all other electrons in the 
system. To shed more light on this point we note that 
at small frequencies uj <^ eV the factor N^^ = eV/\uj\ 
in Eq.® can be associated with non-equilibrium pho- 
ton distribution. This makes the i^coU similar to the 
photocount statistics studied recently by Kindermann et 
al. in Ref. [l9j- In the latter case each photon being 
transmitted through the waveguide and absorbed by the 
photo-detector produces a single count. In the given case 
a single absorption of a photon by the electron gas gen- 
erates the random current pulse in the circuit with a 
zero mean value. The generating function of the current 
distribution in this pulse is given by Il{x,z). Due to 
photon bunching, N^ ^ 1, the pulse is strongly ampli- 
fied thereby producing the long exponential tails in the 
probability distribution P{I). 

"Hot electron" regime, i?Th "^ ^/te- In this linfit the 
collision term in the kinetic equations dominates. There- 
fore the limiting saddle point of the action @ should nul- 
lify the collision integral. To find such solution we note 
that the collision term in the action is invariant under the 
gauge transformation Ge(r) = e~-^«('"^G(;(r)e^"'''-'. Here 
K^{r) = if3[7(r) + /3(r)(e - '/'(r))] and 7,/? and (j) are 
arbitrary functions in space. This leads to the conserva- 
tion of a current density, and a density of the energy flow. 
As well known, the physical Green function G(e, r) with a 
local Fermi distribution f,{r) = [e(^-'l'(''))/nr) ^ij-i nul- 
lifies the collision term in the conventional kinetic equa- 
tion. Its gauge transform Ge{r) does the same for the 
generalized matrix kinetic equation. 

The unknown functions (f), 7, T and P can be found 
from the extremum of the action Fi^ot which is obtained 
by substitution of Ge(r) into the diffusive part of the 
action i^cff- The action i'hot reads 

i^hot = {27ry'gto J dz{-r(V7 - (3V(b)^ (12) 

2 2 

+ (V7 - /3V^)V0 - y t3(V/3)2 + !L(vr2)V/3} 

Here T{z) and (j^iz) are a local temperature and a chem- 
ical potential, while (3{z) and 7(2;) are their quantum 
counterparts. The action H12|l implies boundary condi- 
tions: 0(2)1^^0,= ±e1//2, T{z)\^^^= T, ^7(0) = x and 
7(1) = /3(0) = ^(1) = 0. With the use of integrals of 
motion the Lagrange equations of this action can be re- 
duced to two coupled second order differential equations 
for T{z) and f3{z). We are not aware of their analyti- 
cal solutions under non-zero x and solved them numer- 
ically. The results for the probability distribution P{I) 
are shown in Fig. 2. As in the previous section we have 



evaluated it with the use of the saddle point method. 
Fig. 2 shows that the probability of positive current fluc- 
tuations, A/ > 0, is enhanced in the "hot electron" limit 
as compared to the coherent regime, while the probability 
of negative fluctuations, A/ < 0, is affected in the lesser 
extent. We also note that the action H12|l is equivalent 
to the actions of Ref. |l2L |l3j , which were derived with 
the use of Boltzmann-Langevin approach. These actions 
transforms to i^hot under appropriate change of variables. 

To conclude we investigated the effect of Coulomb in- 
teraction onto the FCS in the one- and two-dimensional 
diffusive conductors. We have revealed the long exponen- 
tial tails in the probability of non-equilibrium big current 
fluctuations in short conductors with i?Th 3> I/t^, te 
being the energy relaxation time. These tails arise from 
the huge fluctuations of the current of electron- hole pairs 
which arc excited by the low frequency fluctuations of 
the electromagnetic field, produced by all other electrons 
in the system. 
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